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We have extended our model for charged current neutrino-nucleus interactions developed in 
Phys. Rev. C 73, 065502 (2006) to neutral current reactions. For the elementary neutrino-nucleon 
interaction, we take into account quasielastic scattering, A excitation and the excitation of the 
resonances in the second resonance region. Our model for the neutrino-nucleus collisions includes 
in-medium effects such as Fermi motion, Pauli blocking, nuclear binding, and final-state interactions. 
They are implemented by means of the Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) coupled- 
channel transport model. This allows us to study exclusive channels, namely pion production 
and nucleon knockout. We find that final-state interactions modify considerably the distributions 
through rescattering, charge-exchange, and absorption. Side-feeding induced by charge-exchange 
\Q • scattering is important in both cases. In the case of pions, there is a strong absorption associated 

f"*) ' with the in-medium pionless decay modes of the A, while nucleon knockout exhibits a considerable 

enhancement of low energy nucleons due to rescattering. At neutrino energies above 1 GeV, we 
£SJ 1 also obtain that the contribution to nucleon knockout from A excitation is comparable to that from 

^ ' quasielastic scattering. 
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> ■ I. INTRODUCTION 
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Neutrino interactions are classified as charged current (CC) or neutral current (NC) processes depending on whether 
a W or a Z boson is exchanged. In the first case, a charged lepton is emitted, whereas the neutrino preserves its 
nature in the second one. The existence and nature of neutral currents played an important role in the establishment 
\Q , of the Standard Model of electroweak interactions. Unfortunately, the experimental study of NC neutrino interactions 
is a demanding task due to the considerable difficulties of collecting data on reactions with cross sections even smaller 
than those of CC processes, and in which the outgoing neutrino leaves no signal, so that the event identification has 
to rely on the detection of one or more hadrons. Still, measurements of NC quasielastic (QE) scattering and pion 
i— h ' production were performed at BNL [l[ , ANL 0, [I[ and Gargamelle 0, Q using deuterium or heavier targets such as 
carbon, aluminum or a propane-freon mixture. 

Nowadays, the discovery of neutrino oscillations has renewed the interest in a better determination of the neutrino- 
nucleus cross sections, aimed to achieve a better understanding of neutrino fluxes and background processes at current 
and future experiments. In particular, it is well known that NC 7r° production in the forward direction is a relevant 
source of background for v e appearance experiments. In fact, both K2K [6] and MiniBooNE 7} have recently measured 
this reaction and the proposed MINERVA experiment plans to do it in the future [8(] . At E v < 2 GeV, pion production 
proceeds mainly through resonance excitation, predominantly of the A(1232) resonance [1, [hJ. But when pions are 
produced in the nucleus, their final-state interaction (FSI) with the nucleons (elastic and charge-exchange scattering, 
absorption) is an essential ingredient of any realistic theoretical description of this reaction. In Ref. [llj . NC pion 
production in nuclei was investigated. In that model, pions, initially produced via A(1232), Pn(1440), and S'ii(1532) 
excitation, undergo a random walk through the nucleus where the pions can change directions but not energy. 

NC neutrino-nucleus interactions are also relevant to answer a fundamental question of hadronic structure, namely, 
the strange-quark contribution to the nucleon spin. Purely isovector CC processes do not depend on the strange 
form factors. Since a non-zero strange axial form factor changes the NC QE cross section on protons and neutrons 
in different ways (see Sec. Ill Al for more details), the ratio R(p/n) of these two cross sections is very sensitive to the 
strange spin, as pointed out by Garvey et al. [12J . Due to the technical difficulties of neutron detection, the proposed 
FINeSSE experiment plans to measure the neutral to charged current ratio i?(NC/CC). In any case, the study of those 
ratios involves both neutrons and protons so it must be performed using nuclear targets. This fact has stimulated a 
considerable amount of theoretical work aiming at the description of nuclear effects in NC nucleon knockout reactions. 
The relativistic Fermi gas description of the nucleus has been adopted in many calculations [HI, [3, [H, [r| ; others 
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use wave functions for the bound nucleons, obtained in relativistic mean- field models [13, EM 0, H3|- The shell 
model [TElfl^l and the continuum random-phase approximation [21, I22I have also been applied. The effect of mesons 
exchange currents has been evaluated by Umino and collaborators 23[. Furthermore, the input from scaling analysis 
of electron scattering data has been used to obtain NC cross sections [24]]. Final-state interactions of the knocked 
out nucleons is neglected in several studies [H, 0, Gil while many treat them with the distorted wave impulse 
approximation [l5lTl^[l^[l8L[20l| and with a multiple scattering Glauber model at higher energies [isj |. The problem 
is that in those approaches it is not possible to take into account nucleon rescattering leading to energy losses, charge 
exchange and multiple nucleon emissions. This can be achieved in Monte Carlo models [25j . It is widely accepted 
that nuclear effects cancel for the ratios of cross sections. However, this is not the case for side-feeding effects caused 
by charge exchange scattering, which are not negligible if the elementary cross section on protons and neutrons differ. 
Such changes in the ratios caused by FSI were obtained by Nieves et al. [25| and we will confirm them here. All 
these theoretical models consider only QE processes, but one should bear in mind that for neutrinos of ~ 1 GeV, 
the excitation of resonances also contributes to nucleon knockout; those resonance excitation events where a nucleon 
is emitted but there are no pions in the final state (or they are produced but not detected) represent a source of 
background for the strange axial form factor measurements, which should be well understood. 

In a recent article (26| we have studied CC neutrino-nucleus interactions in the region of the QE and A peaks. 
Here we extend the model to the NC case. There are three main ingredients in our model: elementary processes, in- 
medium modifications and FSI. QE scattering and A excitation on the nucleon are treated in a relativistic formalism, 
incorporating state-of-the-art parameterizations of the form factors for both the nucleon and the N — A transition. It 
is, however, known that the excitation of heavier resonances is not negligible above (» 1,5 GeV d, (hJ. Therefore, 
in this article, we also consider the excitation of the N* states Pn(1440), D 13 (1520) and S , n(1535), by means of the 
Rein and Sehgal model [T3| , which has been extensively used in the simulations of neutrino experiments [13, HH, • 
We should, however, recall that new information on the N — N* electromagnetic transition form factors is available 
from the analysis [30( of recent electron scattering data. Hence, there is room for improvement in this part of the 
model [31,]. Next, we take into account nuclear effects: Fermi motion, Pauli blocking and the binding of the nucleons in 
a density and momentum dependent mean-field. Finally, FSI are implemented in the framework of the semiclassical 
coupled-channel transport theory: the Giessen BUU model. With these tools we investigate pion production and 
nucleon knockout in neutral current neutrino-nucleus interactions at intermediate energies. 

Next, we present our model for the elementary cross sections emphasizing the structure of the neutral hadronic 
currents and the role of strange form factors. Following an overview of the nuclear model and final-state interactions, 
we present our results and compare these to other calculations. Summary and conclusions are given at the end. 



II. NEUTRINO-NUCLEON REACTIONS 



In this section we present the model we have adopted for the description of elementary neutrino-nucleon interactions, 
emphasizing the aspects that are specific of neutral current processes. Further details on our corresponding approach 
to charged current processes can be found in Ref. [26| . 

We consider neutral current reactions of the type, 

v{k) + N(p)^v(k')+X(p'), (1) 

with k a = \ E U) k\ k' a — (e' v , fc'^j, p a — (E, p) and p' a — (e', p'j. The cross section can be cast as an integral 
over the azimuthal angle of the outgoing neutrino 



dV„jv f J , 1 1 1 



dQ 2 d£^ J 64tt 2 \k-p\ E v 



p l2 -M' 2 )\M\\ (2) 



where Q 2 — — (k — k') 2 and M' is the (invariant) mass of the outgoing baryon. For neutral current (NC) interactions, 
the matrix element squared, summed and averaged over spins, 

\M\ 2 = ^fL a pH^, (3) 

is given in terms of the Fermi constant Gf = 1.16637 x 10~ 5 GeV -2 , the leptonic tensor L a p and the hadronic tensor 
H a @ . While the calculation of L a p is straightforward, the hadronic current entering H a P has to be parametrized in 
terms of form factors and thus depends on the specific reaction. Two processes dominate the reaction at neutrino 
energies up to about 1.5 GeV, namely quasielastic scattering, 

vn — > vn, vp — > vp, (4) 
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and A(1232)P 33 excitation, 

z/n— >i/A°, i/p— >z/A+. (5) 

Here, in addition, we consider the excitation of the N* states R = Pn(1440), £>i3(1520) and Sii(1535), which form 
the so-called second resonance region, 



vn 



vR°, vp^vR+. (6) 



A. Quasielastic scattering 

The hadronic tensor H a P for quasielastic scattering is determined by the hadronic current J® E given by 

J$ E = (N'\J» c (0)\N) 

= (N'\(V a NC -A^ c )(0)\N) 

= u{p')B a u{p) (7) 

with q a — p' a — p a , N = p,n and 

B a = ( 7a - *p) F? + ^^FT + 7a rf + fi^F?, (8) 

where Mjy denotes the nucleon mass; F^ 2 , Fa an d Fp are * ne vec tor, axial and pseudoscalar form factors, respectively. 
Due to time invariance, these form factors are real functions of Q 2 . In the expression for the cross section, Fp 
appears only multiplied by the neutrino mass so we ignore it from now on. The term {^q a )/q 2 ensures vector current 
conservation even if the masses of the initial and final nucleons differ, as might be the case in the presence of a 
momentum-dependent nuclear mean-field potential. The hadronic tensor H a P then follows as 



H oe = ^ Tr 

QE 2 



(fl + M)B a (f +M')BP (9) 



with 

Ba = loBl 70 (10) 

and — in the free nucleon case with M = M' — M^. 
The vector current has the form 

V a NC = (1 - 2sin 2 M^ 3 -2s^e w \jl lj' a , (11) 

where Q\y is the weak-mixing angle (sin 2 9w = 0.2228), is the third component of the isovector current, is the 
isoscalar (hypercharge) one and J a stands for the strange part. All of these terms have the same Dirac structure. 
Therefore, the NC form factor can be written as 

K 2 = Q-2sin 2 ^Ff >2 -i^ 2 -l^ 2 , (12) 

Fh = (I ~ 2 sin2 <V) F h l F i,2 ^1,2 (13) 
in terms of the the standard Dirac and Pauli form factors of the nucleon Ff' n and F 2 ' n and a strange component 
Analogously, the axial current consists of the third component of the isovector axial current and a strangeness part 

A™ = Al + \A> a . (14) 



This implies that 



F% n = ±-F A + -F%, (15) 
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where Fa is the axial form factor for charged current QE scattering. Notice that here it appears with different signs 
for protons and neutrons. Fj± stands for the strange axial form factor. 

As in our earlier work f2(|, we use the BBA-2003 parametrization [32I ] for the non-strange vector form factors and a 
dipole ansatz with Ma = 1.0 GeV for Fa (See Ref. 33:] for an overview on the extraction of Ma from neutrino-nucleon 
scattering, pion electroproduction and muon capture). 

The strangeness content of the nucleon, encoded in Ff 2 and is still an open question. It can be investigated 
in a combined study of parity- violating polarized electron scattering and neutral current neutrino scattering. Parity- 
violating electron scattering is very sensitive to the strange electric and magnetic form factors (i. e., to the strange 
vector form factors) and much less to the strange axial vector form factor (cf. e. g. Ref. (34|). The opposite holds 
for NC neutrino scatterin g. A n extensive program on parity violation has evolved in the last years: the SAMPLE 
experiment at MIT /Bates |35j. HAPPEX (3(|[37[ and GO [38| at JLab, and PVA4 ^ in Mainz have extracted linear 
combinations of the strange electric and magnetic form factors at different Q 2 values. Recently, the strange electric 
and magnetic form factors were extracted from a combined set of available parity-violating electron scattering data 
by Young et al. [1(3]. However, data on NC neutrino scattering, needed to determine the strange axial form factor, 
are scarce. The best measurement to date is the E734 experiment at BNL It measured neutrino-proton and 
antineutrino-proton elastic scattering albeit with large systematical errors and only small statistics. Former attempts 
to extract the strange axial form factor from the data jU |42j faced the fact that, as pointed out by Alberico et al. f42| . 
"the experimental uncertainty is still too large to be conclusive about specific values of the strange form factors of 
the nucleon" and a "rather wide range of values for the strange parameters is compatible with the BNL E734 data." 
The advent of new polarized electron scatte ring data from the above mentioned experiments changes appreciably 
the situation because, as shown by Pate [43|, |44j . it allows to perform a simultaneous determination of all (electric, 
magnetic and axial) strange form factors with small error bars, in spite of the uncertainties of the E734 data [45|]. But 
this is only possible in the region of 0.45 < Q 2 < 1.05 GeV 2 where the E734 differential cross sections were measured, 
so new NC neutrino scattering data at low Q 2 are needed for a reliable extrapolation down to Q 2 = 0. The proposed 
FINeSSE experiment [46| will hopefully fill this gap. 

In view of the low sensitivity of NC neutrino scattering to the strange vector form factors and the large sensitivity 
to the axial one (see Ref. [47J for a detailed study using the parameterizations of Garvey et al. [41() here, for the sake 
of simplicity, we choose to take 



and 



Ff(0) 

nm = 



= 0, 
= 0, 

As 



1 



Mi 



(16) 
(17) 

(18) 



assuming that the strange axial mass is equal to the non-strange one. Here As denotes the strange contribution to 
the nucleon spin. In line with Ref. (l5| . we use As = —0.15 and As = as representative values. 

The cross sections for NC QE scattering on proton and neutron are shown in Fig. [TJ Note how the strange spin 
causes opposite effects on the cross sections for protons and neutrons; this is a direct consequence of the two different 
signs in Eq. (fl"5|) . 



B. Production of the A and of higher resonances 

The hadronic tensor H a ^ for A excitation is determined by the hadronic current given by (46 



with 



J. 
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(19) 
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(20) 
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FIG. 1: Differential and integrated cross section for NC quasielastic scattering on protons (top) and neutrons (bottom). The 
solid lines denote the results with As = —0.15, the dashed lines denote the ones without the strange axial form factor (As = 0). 



where ip (p') is the Rarita-Schwinger spinor for the A, and u(p) is the Dirac spinor for the nucleoli. The hadronic 
tensor follows to 



H 



a/3 _ 



i T r 



M)B ap A pa B aP 



with 



B a f3 = JnBl^o 



and, for free nucleons, M = Mj\f. The spin 3/2 projection operator is given by 



A„ 



9pa 



2 p'p: 



pro 



1 V'plo ~ Palp 1 



7pJ<T 



(21) 
(22) 

(23) 



In contrast to quasielastic scattering, where the neutral current is sensitive to the isoscalar quark content of the 
nucleon, the TV — A transition is purely isovector. Therefore, the neutral current reduces to 

J^ c = (l-2 S m 2 9 w )V*-Al, (24) 

where and A„ are the third components of the vector and axial isovector currents, respectively. Thus, we have 



Cf = (1 -2 sin 2 6W)C, 
Cf = Cf, 



v 

i ' 



(25) 
(26) 



where C^and C? are the vector and axial charged current transition form factors for which we use the parametrization 
of Ref . I49| as in our earlier work [291 . 
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FIG. 2: Differential and integrated cross section for NC resonance excitation on protons (top) and neutrons (bottom). The 
cross section for A excitation (solid lines) clearly dominates. The lines denoting the cross section for D13 overlap with the ones 
for the S11. 



The phenomenological information about the neutrino induced N — N* transition is far more scarce. Several articles 
have considered N — N* vector form factors derived from helicity amplitudes @, [E(J, [EH, [HJ , extracted from electron 
scattering experiments or calculated with different quark models (ioj . especially for the Pn(1440) [53l. IHiL IHR l56l . 
In the context of neutrino scattering, the most recent study based on new JLAB electron scattering data has been 
performed by Lalakulich et al. [31]. Experimental information on the N — N* axial form factors is very limited. 
Goldberger-Treiman relations have been derived for the axial couplings [3, [E3] , but there is no information about the 
Q 2 dependence. For the production of higher resonances, we use the matrix elements derived from the model of Rein 
and Sehgal [lfj who apply a quark model to calculate the vector and axial N — N* transitions. 1 This model has been 
widely used in the simulation and analysis of many neutrino experiments [27L |28[ . |29| . The Q 2 dependence of their 
form factors has the form of a modified dipole (see their Eq. (3.12)). For the dipole masses, we use Ma — 1.032 GeV 
and My = 0.84 GeV (cf. e. g. Ref. [58|]). In this way we include the production of the resonances Pn(1440), Z?i3(1520) 
and S n (1535). 

The finite width of the resonances is accounted for in the cross section of Eq. ([2]) by replacing 

5 (V 2 - M' 2 ) - -lira (- i =-) , (27) 

where M' is the pole mass of the resonance and T is the energy dependent full decay width in the vacuum. We use 
a consistent set of resonance parameters taken from a single analysis, namely the one of Manley and Saleski (59|. All 
the relevant decay channels as ttN, it A, pN and r]N are included. 



1 Our matrix element M corresponds to T in the notation of Rein and Sehgal fioll . 
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FIG. 3: Integrated cross section for NC one-pion production via resonance excitation and subsequent decay. The left panel 
shows the pion production on protons and the right one on neutrons. Data are shown for the reaction vn —* vpn~ (taken from 
Ref. 0]). 
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The cross sections for NC resonance excitation on protons and neutrons are shown in Fig. [5] In both cases, the A 
yield is considerably larger than the other ones. 

In the vacuum, we model pion production via resonance excitation and their subsequent decay (we emphasize that 
in the nuclear medium resonances can undergo "pionless decay" (cf. Sec. IIII A")) ). When a neutrino scatters with a 
nucleon, either an isospin 3/2 resonance (A) or an isospin 1/2 resonance {P\\, -D13 and Su) is produced. They decay, 
among other channels, with a certain branching ratio into irN pairs. By including the appropriate Clebsch-Gordan 
coefficients we obtain for the cross sections of the four possible one-pion production channels 

12/ \ 

3^a+ + 3 (M P + + b 2 a D + + b 3 (T s +J , (28) 
2 1/ \ 

3°A+ + 3 (^i°>+ + &2<7 D+ + b^s+J ' ( 29 ) 
12/ \ 

3^ao + 3 (bivp^ + b 2 <r D o 3 + has^J , (30) 
2 1/ \ 

3^A« + 3 [bicrp^ + b 2 (T D a 3 + b 3 a s oJ , (31) 

where (ctro) are the cross sections for resonance excitation on protons (neutrons) as shown in Fig. [2] For 
consistency, the branching ratios bi are also taken from the analysis of Manley and Saleski [59| which gives b\ = 0.69, 
b 2 = 0.59 and 63 = 0.51. In Fig. [3l we plot the results for the integrated pion production cross section on protons (left) 
and neutrons (right). Data on NC pion production is extremely sparse. A measurement was performed on D2 at the 
ANL bubble chamber for the vn — > vpn~ channel @ — these data are also shown in Fig. f3J The remaining channels 
have only been measured at the Gargamelle bubble chamber Q (see also the reanalysis by Hawker 60] ) on a propane- 
freon mixture and not on "elementary targets". Thus, we do not show them here. However, we have compared our 
results to the calculations of Paschos et al. [ll[ (cf. their Figs. 4-7) and have found reasonable agreement for the 
integrated vacuum one-pion production cross sections. 



III. NEUTRINO-NUCLEUS REACTIONS 
A. In-medium modifications and final-state interactions 

Nuclear effects play the central role in this study. A detailed description of our nuclear model is given in Ref. [26| 
where we applied it to CC neutrino-nucleus scattering. Here we only outline the main features. 

The cross section for neutrino scattering on free nucleons, given in Eq. ([2]), has to be modified when the nucleon is 
bound in the nucleus. First, we consider Fermi motion of the initial nucleons and Pauli blocking of the final ones in 
a local density approximation, in which the local Fermi momenta of the nucleons are given by 

/3 \ 1/3 

p F {r)=[--K 2 p{r)\ . (32) 
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For the density distribution of heavier nuclei we use a Woods-Saxon form with parameters extracted from Hartree-Fock 
calculations. For 12 C, we take a harmonic oscillator density as given in Ref. [6l|. 

The nucleons are bound in a density and momentum dependent scalar potential U(r, p), whose parametrization 
has been obtained by fits to the saturation density of nuclear matter, and also to the momentum dependence of the 
nucleon optical potential as measured in pA collisions [62j • Since the A is less bound in the nucleus than the nucleons 
- at normal nuclear density one finds Ua ~ —30 MeV versus Un ~ —45 MeV [63| — , we set the A potential to 
2/3 of the nucleon potential. For the other resonances the same potential as for the nucleons is adopted. A particle 
bound in the nucleus acquires an effective mass, defined as 

M eff = M N>R + U NjR { f , p ), (33) 

where r is the corresponding "free" mass. As consequence of the density and momentum dependence of the scalar 
potential the effective masses of initial and final particles in a scattering process can be different even if their rest 
masses are equal. We account for this fact by replacing M and M' (note that Mjv remains unchanged) in the cross 
section formula (Eq. ((2|)) and in the hadronic tensors (Eq. ([9]) and Eq. (|2Tj) ) by their respective effective masses (cf. 
Ref. [H for details). 

The final-state interactions (FSI) of the produced particles are implemented by means of the coupled-channel semi- 
classical Giessen Boltzmann-Uehling-Uhlenbeck (GiBUU) transport model [64]. Originally developed to describe 
heavy- ion collision s l65l. 66 ,l67l , it has been extended to describe the interactions of pions, photons and electrons with 



nuclei H M M MMW- 

Recently, we further extended the GiBUU model to describe neutrino scattering off 

nuclei [26j. This extension does not require the introduction of any new nuclear parameter. 

In our model, the space time evolution of a many-particle system under the influence of a mean-field potential 
and a collision term is described by a BUU equation for each particle species. A collision term accounts for changes 
(gain and loss) in the phase space density due to elastic and inelastic collisions between the particles and also due to 
particle decays into other hadrons whenever it is allowed by Pauli blocking. In between the collisions, all particles (also 
resonances) are propagated in their mean-field potential according to their BUU equation (cf. Refs. [H,[7(|). More 
details on the cross sections, their in-medium modification and the included particles are given in the aforementioned 
references. 

Inside the nuclear medium, the width of the resonances is modified. Since their decay into ttN pairs might be Pauli 
blocked, their width is lowered compared to the vacuum width. On the contrary, we have an increase of the width 
due to collisions in the medium. Elastic and inelastic scattering of the resonances with the nucleons contribute to this 
collisional broadening. Therefore, the total in-medium width is given by a sum of the Pauli modified decay width T 
and the collisional width V co u- Despite the decrease caused by Pauli blocking, the dominant effect is a broadening of 
the resonances in the medium. 

The most relevant states for neutrino-induced reactions at intermediate energies are the nucleon, the A resonance, 
the pion and their interactions. For the NN cross section and its angular dependence we use a fit to data from 
Ref. [73| . For the pion cross sections we use a resonance model with the background fitted to data as shown in detail 
in Ref. [7l|. The A resonance is propagated off-shell in our model and decays isotropically. Its decay into a pion 
nucleon pair is Pauli blocked if the momentum of the nucleon is below the Fermi momentum. We allow not only for the 
decay of the A, but also for FSI of the A in the nuclear medium. Possible absorption processes include AN — > AW, 
ANN NNN and AN — > ttNN. All those processes are implemented by an absorption probability depending on 
the in-medium collisional width for which we use the results of Oset and Salcedo (74J. We emphasize that the whole 
ttN A dynamics has been tested extensively and also compared to data in A A — > it X [65| , it A — > X [7l|, [75| and 
7 A — > 7r X [6^| and recently in pionic double charge exchange reactions [72| . 

In conclusion, FSI lead to absorption, charge exchange, a redistribution of energy and momentum as well as to 
the production of new particles. In our coupled-channel treatment of the FSI — in which the BUU equations are 
coupled through the collision term and, with less strength, also through the potentials — our model differs from 
standard Glauber approaches because the collision term allows not only for absorption but also for side-feeding and 
rescattering. In addition, our model is applicable in many different nuclear reactions using the same physics input. 
In particular, an important prerequisite for any model aiming at the description of the interaction of neutrinos with 
nuclei is a good description of electron- or photon-induced reactions. There, extensive tests against existing data are 
possible and have been successfully performed (7f| [77J ■ 



B. Results and discussion 



We now present our results for NC pion production and nucleon knockout for neutrino energies up to 2 GeV for 
nuclei commonly used in neutrino experiments. 
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FIG. 4: Integrated cross section for 7r° (top), ir + (bottom left) and 7r~ (bottom right) production on 56 Fe as a function of E v . 
The dashed lines show the results without FSI (only the decay of resonances is possible); the results denoted by the solid lines 
include FSI. Also indicated is the origin of the pions (QE, A excitation or higher resonances (hiRES)). 



1. Pion production 

We begin our discussion on pion production with the total production cross sections for 7r + , tt° and ir~ on 56 Fe which 
are shown in Fig. [3J The dashed line denotes the pions stemming from the decay of the initially produced resonances 
— no further FSI are taken into account here. The cross section for tt° production is significantly higher than those 
of the charged channels. This is a direct consequence of the isospin structure of the resonance decay. "Switching on" 
FSI allows the resonances to interact in different ways besides a simple decay. Also the produced pions can interact 
further or undergo absorption. These FSI lead to a strong reduction of the total yield in the ir° channel (compare solid 
and dashed line in the upper panel). The reduction is much smaller in the ir + and ir~ channel because the 7r° undergo 
charge exchange scattering and thus contribute to the charged channels (side-feeding). The effect in the opposite 
direction is less important due to the smaller elementary tt + and 7r~ production cross section. Thus, side-feeding 
shifts strength always from the dominant into the less dominant channel. In the context of neutrino reactions, this 
effect was first described by Adler et al. 78] . It is also observed in charged currents from the ir + into the ir° channel 

[2ilzlils3. 

Pions can also emerge from the initial QE neutrino-nucleon reaction when the produced nucleon rescatters producing 
a A or directly a pion. However, as can be seen from Fig. [4] this process is not very sizable because it is relevant only 
at high Q 2 (dash-dotted lines). 

A more detailed understanding of pion production is obtained by studying the pion kinetic energy distributions. 
In Fig. [5] we show the kinetic energy spectra for 7r°, ir + and ir~ production for three neutrino energies. The dashed 
lines denote again the result without final-state interactions and the solid lines the result of the full calculation. The 
contributions from initial A excitation (dotted lines) and from initial QE events (dash-dotted lines) are also shown. 
Pion production through FSI of initial QE processes contributes mostly to the low energy region of the pion spectra 
due to the redistribution of the energy in the collisions. While the overall shape of the dashed lines (without FSI) is 
dictated by the predominant p-wave production mechanism through the A resonance, the shape of the solid lines (full 
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FIG. 5: Kinetic energy differential cross section for tt production on Fe versus the pion kinetic energy T w at different values of 
E v . The dashed lines denote the calculation without FSI where only the decay of resonances is included; the solid lines denote 
the one with FSI. Also indicated is whether the pion comes from initial QE or A excitation (dash-dotted or dotted lines). 



calculation) is influenced by the energy dependence of the pion absorption and rescattering. The main absorption 
mechanism for pions above T„ w 0.1 GeV is irN — > A followed by AN — > NN which leads to a considerable reduction 
of the cross section. Elastic scattering tt N —* irN redistributes the kinetic energies and thus also shifts the spectrum 
to lower energies. In the case of the smaller ir + and ir~ channels the already discussed side-feeding enhances the peak 
in the middle and bottom panels of Fig. [5] over the value obtained without FSI. Note that these spectra are very 
similar in shape to those measured in (7, n°) reactions on nuclei (cf. Figs. 13 and 14 in Ref. [8l|). 

The impact of FSI is even more visible in the ratios obtained by dividing the differential cross section with FSI 
by the one without FSI. These are plotted in Fig. [6] for 56 Fe and 16 O at E v = 1 GeV versus the pion kinetic energy 
(short-dashed for ir + , solid for ir° and dashed lines for 7r _ ). As expected, the absorption is bigger in the heavier 
nucleus ( 56 Fe) than in the lighter one ( 16 0) — it scales with A 2 ^ 3 . For pions with kinetic energy > 0.1 GeV we 
find strong effects of FSI. This is the region where pion absorption and rescattering are most prominent due to the 
excitation of the A resonance around its peak position. At lower energies (TV ~ 0.07 GeV) we find a peak because 
pions of higher energy in average lose energy via rescattering. At still lower pion energies, the multi-nucleon pion 
absorption mechanism takes over, leading to a small dip. We stress that a similar pattern has been experimentally 
observed by Krusche et al. [8l| in pion photoproduction (cf. Fig. 16 in Ref. 81]). This particular dependence of the 
ratio reflects well-known features of the 7riVA dynamics in nuclei. 

It would be interesting to compare these results with earlier calculations of neutrino-induced pion production on 
nuclei by Paschos et al. [Til ] . However, these authors have very recently found an error in their elementary pion- 
production spectra [HJ which affects their Figs. 8-14 for the pion energy distribution in Ref. [ll[. In view of this we 
abstain from a comparison until the corrected results are published. 
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2. Nucleon knockout 



We now continue our discussion on exclusive channels with nucleon knockout. We consider all nucleons leaving 
the nucleus in the vA reaction which — in case FSI are included — is not necessarily a single- nucleon knockout. In 
Fig. [7] we show the integrated cross sections for proton and neutron knockout on 56 Fe. The solid lines, showing the 
results with FSI included, lie in both cases clearly above the ones without FSI (dashed lines). This enhancement 
is entirely caused by secondary interactions. Furthermore, the initial process leading to a knocked out nucleon is 
indicated. In the vN collision either QE scattering (dash-dotted), A (dotted) or N* excitation (double-dashed lines) 
is possible. While the pion cross section was dominated by the initial A excitation, here initial QE and initial A 
excitation contribute to the total cross section above E v ss 1 GeV with nearly equal amounts. This reflects the energy 
dependence of the vN cross sections (cf. Fig. [T]and[2]). Only up to neutrino energies of w 0.5 GeV the resonance 
contributions to nucleon knockout is negligible. We note that our method of using the A selfenergies of Oset and 
Salcedo [7|[ (cf. Sec. IIII A|) even underestimates the number of knocked out nucleons because we do not treat the 
process A7V — * NN cxplicitely. Our A resonance contribution of knocked out nucleons stems solely from the decay 
A — > ttN . 

In Fig. [8] we plot the kinetic energy differential cross section for proton and neutron knockout versus the kinetic 
energy on 56 Fe for different values of E u . The inclusion of FSI strongly modifies the shape of the distribution (compare 
the dashed and the solid lines). This effect is due to the rescattering of high energy nucleons in the medium which 
reduces the flux at higher energies while, simultaneously, a large number of secondary nucleons at lower energies is 
emitted. At E u = 0.5 GeV, nucleon knockout is clearly dominated by QE processes, however, the plots show how the 
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FIG. 8: Kinetic energy differential cross section for proton (left) and neutron (right) knockout on 56 Fe versus the nucleon kinetic 
energy Tjv at different values of E v . The dashed lines denote the calculation without FSI where only the decay of resonances 
is included; the results denoted by the solid lines are with FSI. Also indicated is whether the pion comes from initial QE or A 
excitation (dash-dotted or dotted line). 



A becomes progressively more important as E v increases. 

In contrast to the CC reaction studied in Ref. [2(|, where neutrons were only emitted by secondary collisions, in 
the NC case, both the neutron and proton kinetic energy distributions are equally affected by FSI since their total 
yields without FSI are comparable. This can be seen also in Fig. [9] where we plot the ratio of the proton to neutron 
kinetic energy differential cross section for 56 Fe and 12 C at E u = 0.5 GeV. The calculations with and without FSI 
(solid and dashed lines) agree approximately if As = —0.15. This shows that the effect of the final-state interaction 
cancels out regardless of the nucleus. This is different in the CC case (cf. Ref. 26]) where we found strong effects 
of sidc-fceding from the dominant p channel into the suppressed n channel. In fact, side-feeding is only important 
when the initial proton and neutron yields are different. Indeed, when we "switch off" the strange axial form factor 
(i. e. take As — 0) then the elementary proton and neutron yields are different as shown in Fig. [TJ while the neutron 
cross section is enhanced, the proton cross section is reduced. Therefore, we expect side-feeding from the neutron to 
the proton channel which changes the p/n ratio as observed in Fig. [5] (dash-dotted versus dotted lines). Since nucleons 
lose energy in (charge exchange) scattering, the effect is more pronounced at low kinetic energies. 
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Finally, we compare our results on nucleon knockout with other calculations. In Fig. [10] we show our results for 
the kinetic energy differential cross section together with those of Nieves et al. [25| (denoted with "NVV"). Since 
they do not include any resonances, we have "switched off" the initial resonance excitation in our calculation, so only 
nucleon knockout induced by initial QE events is considered. The discrepancy of our result without FSI to the result 
without NN rescattering of Ref. [25[ (dashed versus dotted line in Fig. [10]) could be attributed to — in addition to 
differences in the momentum distribution and the potentials — the RPA correlations included in their calculation 
which lead to a reduction of the cross section and a spreading of the spectrum. To model the rescattering of the 
primary nucleons in the nucleus, Nieves et al. use a Monte Carlo simulation with elastic NN cross section similar to 
ours. 2 Therefore we expect a similar behavior when FSI are included. Indeed, as one can see, when the rescattering of 
the outgoing nucleons is "turned on" , both calculations lead to very similar results, namely a reduction of the flux for 
higher energetic nucleons and a large number of secondary low energy nucleons (in Fig. [TU] solid versus dash-dotted 
line). Also for the p/n ratios we find reasonable agreement: our ratios plotted in Fig. |2]show a behavior similar to 
the ones of Nieves et al. (cf. right panel of Fig. 17 in Refs. [HI). 

Furthermore, we compare to the calculation of Martinez et al. 18]. While our results without FSI agree approx- 
imately with their relativistic plane wave impulse approximation (RPWIA) (up to small differences which can be 



2 We emphasize that, in addition, we allow for inelastic NN collisions. 
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attributed to the momentum distribution of the nucleons (relativistic mean-field versus local Fermi gas) and to the 
potentials used in the calculations), a comparison of the results with FSI may be meaningless, because our method 
describes exclusive nucleon knockout reactions which may contain more than one nucleon in the outgoing channel and 
may leave the residual nucleus in a highly excited state. On the contrary, in their model, the FSI of the nucleon are 
considered within two frameworks. In one case, they use a relativistic optical potential in a distorted-wave impulse 
approximation (RDWIA), in the other case a relativistic multiple-scattering Glauber approximation (RMSGA) is 
used. These models can therefore account only for single-nucleon knockout. We emphasize that these absorption 
mechanisms as used by Martinez et al. can describe only the flux reduction at higher nucleon energies but not the 
rescattering in the medium which leads to the emission of a large number of lower energy secondary nucleons. Nu- 
cleons are not just absorbed but — through rescattering — ejected with a different energy, angle and/or charge. In 
addition, we stress that already at E v ks 1 GeV a large par t of the ejected nucleons stems from A excitation (see 
Fig. [7]), a process not contained in the calculation of Rcf. 18]. 



IV. SUMMARY AND CONCLUSIONS 



In the present work we have investigated NC neutrino interactions with nucleons and nuclei extending our earlier 
work on CC reactions (26[. The neutrino- nucleon reaction is dominated by quasielastic scattering and A excitation; 
we have found that the higher resonances only give a minor contribution at neutrino energies below 1.5 GeV. 

Nuclear effects are included in the framework of a coupled-channel BUU transport theory (GiBUU model) where we 
account for in-medium modifications from Fermi motion, Pauli blocking, nuclear binding and collisional broadening 
of resonances. In this model, FSI are implemented by means of transport theory which permits exclusive reactions 
to be studied. Within the same model, we can describe neutrino-induced pion production and nucleon knockout. In 
this respect, the model presented here is unique. 

The pion production cross section is especially influenced by final-state interactions. In the elementary neutrino- 
nucleon reaction, more ir° than charged pions are produced due to the isospin structure of resonance decay. When 
final-state interactions are included, those tt° get absorbed or reinteract leading to side-feeding in the smaller tt^ 
channels. Quasielastic scattering followed by 7r production in NN collisions also accounts for a small fraction of the 
pion production cross section. 

Also for nucleon knockout the influence of the final-state interactions is significant. The rescattering of high energy 
nucleons in the nucleus leads to a reduction of higher kinetic energy nucleons while a large number of secondary 
nucleons at lower nucleon kinetic energies are ejected. Also in the case of nucleon knockout, we found that side- 
feeding is important. Furthermore, we have illustrated that for neutrino energies > 1 GeV initial resonance excitation 
(predominantly A) leads to a significant contribution to nucleon knockout. 

Summarizing, we have found that in-medium modifications, and especially final-state interactions, have a big 
influence on the neutrino-nucleus cross sections. We emphasize that a good and well tested description of these effects 
is crucial for the understanding of current and future experiments. 
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